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ABSTRACT 

The technique of positron amission tomography PET 
promises to open up naw medical frontiers by enabling physicians 
to stxady the metabolic activity of the body in a pictorial 
manner* A suitable radio~active isotope [3 ] / which emits 
positrons with a short half life time, is injected into the 
patients body* The radio-isotops flows with the blood stream 
and distributes itself in the various organs of the body* Tho 
positrons emitted by the radio-isotops interacts with the 
electrons in the nuclei of the surro^anding material and get 
annihilated in the process* This annihilation produces a pair 
of photons which are sxibsequently detected in coincidence by 
an array of discrete detector elements surrounding the body* 

The short physical halfr-lifa of uhe majority of the positron 
emitting radio-nuclide lowers ths radiation dose to the patient 
and permits repeat studies at short time interval [3]» Much 
as in X-Ray transmission tomography and other modes of 
camputerizad tomography ths quality of the racons time ted imago 
is very sensitive to the mathematical algorithm to be used for 
reconstruction* The objective of this study is to compare the 
performance of the EM algorithm [8] (specifically designed 
for PET) and an adaptation of the convolution bach projection 
algorithm [ 7 ] widely used in transmission tomograph CAT scanner* 



Chapter 1 


A discussion of 
Positron emission tomography 

Positron emission tomography (PET) is a medical diagnostic 
technique, which enables a physician to study blood, flow and 
metabolic activity of an organ, in a visual way# For this, a 
biochemical metabolite labelled with a positron emitting 
radioactive material is introduced into the organ, under study. 
The radioactive emissions are then counted by using PET scanner- 
a machine consisting of detector elements mounted on one or 
more rings, positioned so that it surrounds the patient's body. 

The choice of biochemical and the radioactive tracer 
depends on the organ to be studied and on the questions of 
interest# For example, labelled glucose is often used in PET 
studies of the brain while metabolism of the heart is measur'sd 
with labelled deoxyglucose and labelled palmitic acid# 

PET scans of brains of people suffering from schizophrenia 
and certain other diseases shox'/ distinctive metabolic patterns 
connected with these dis';asos# These metabolic portraits 
obtained from PET scans plays an important role in diagnosing 
such diseas =;s and in assessing the effectiveness of various 
treatments# 

An overview of the subject, including a discussion on 
clinical, biochemical, and other interesting aspects of PET 
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can be found in Brownell et al# (1932) and in Ter-Pogossian 
et aU (1930). 

THE physics of THE PET : 

When a positron is emitted from a radioactive positron 
emitting material^ it gets annihilated with a nearby electron* 

The an::ihilation creates two X-Ray photons that fly off the 
point of annihilation, at the spaed of light, in (nearly) 
opposite directions, along a line with a completely random 
orientation* This pair of X-Ray photon is detected in coincidanci 
by a pair of detector elements that define a cylinderical volume 
(to be referred to as a detector tube or, sirrply a tube)* 

The set of data collected in a PET scan is the txibe counts 

n*(l), n*(2),** » ,n*(D), T^fliere n*(d) is the total number of 

"till 

coincidences counted by the d^ detector t\ibe and D is the total 
number of tubes. 

Total tubes count is much smaller than the total ntmber 
of emissions, because 

(1) Some of the photons do not get caught by detector ring* 

(2) Some of the photons are absorbed by the body'^s tissue. ! 

These facts are schematically described in fig. (1*1)* 
Table (1) shows the principal radioisotopes in diagnostic use. 
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CHOICE OF DETECTOR MATERIALS : 

Three msterials are ^andsr considerations for use in 
positron tomographs : 

(1) Sodium Iodide (Nal(Tl)) 

(2) BISMUTH GERRZ^NATE (Bi^ Ge^ 0^^^* (BGO) (Cho and Parukhi 1977 

(3) Caesium Fluoride (CsP) 

(Allemand et al* (l9B0). 

All these three have widely differing properties which are 
summarized in the table (2)* 

BGO is twice as dense as Nal(Tl) which results in a very 
high detection efficiency even for small crystals* CsF has 
excellent coincidence resolving time. Fur-thermore, BGO is not 
hydroscopic and therefore, does not need to be hermetically 
sealed as does Nal(Tl) and CsF. Most high resolution devices 
with detectors below 10 mm in size have BGO detectors (Deranzo 
et al. 1931 J Erikson at al« 1930 J Nohra at al* 1980 | Brooks 
at al# 1980) because of the high detection efficiency and good 
packing fraction* CsF is preferred for lower resolution 
devices intaaded for fast dynamic studies (Yamamoto et al 1932). 
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a 



b 

figure I. (a) Two anrMa^oos: one, at x, that is detected in tube 
P, $7) and the other, at jr, that passes undetected because the photon 
path does not intersect the detector ring, (b) Top view of the detector 
rmg *piane* of a. 

Ref. [9 ] 
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Table 1.-. Principal radioisotopes in diagnostic use (excepting C, N, O, F) 


Isotope 

Radiophannaccut ica! 

Use 

Rcfcrcncsc 

“Ga 

Proteins 

Regional blood-volume 

Mazi&re, B. ei al (1982) 

• 

Microsphcrcs 

Regional blood-liow 

Yvert et al. (1979) 

"Co 

Bleomycin 

Tumour tracer 

Lagunas-Solar and 
Jungerman (1979) 


DTPA 

CSF kinetics 

Mazi^re, B. «f al (1981) 

“Rb 

— 

Myocardial blood- 
flow 

Ltmbre^l and Wolf 
(1979) 

w.76Br 

Bromocriptine 

Dopamine receptors 

Huang ef al, (1981) 


Bromospiroperidol 

Brain studies 

Friedman « al (1981) 

’»Se 

Methionine 

Protein metabolism 

Guillaume et al (1979) 


Chloride 

Myocardial imaging 

Ku 0 al. (1979) 


Chelate 

Kidney function 

Atdicr 0 al (1980) 

I23j 


Cardiac angiography 

Richards and Ku (1979) 

>»Ne 

— 

Local lung ventilation 

Crouzcl 0 a!. (1980a.b) 


Ref. [3, page 21 J 


Table 2 , . Comparison of some detector materials 



NalCri) 

BGO 

OF 

Density (gem”’^) 

3.67 

7.13 

4.61 

Atomic number 

11,53 

82,32,8 

55,9 

Relative light output 

1.0 

0.08 

0.05 

Scintillation peak 

410 

480 

390 

wavelength (nm) 

Index of refraction 

1.78 

2.15 

1.48 

Scintillation decay 

230 

300 

5 

time (ns) 

Time resolution** 

-12 

-20 

-6 

(2t, ns) 

Detection cfficicncyt 

Ciy^al width 

20 mm 

61(41)% 

90(88)% 

68(49)% 

8 mm 

52(31)% 

84(79)% 

57(38)% 

4mm 

46(25)% 

78(71)% 

50(29)% 


• Coincidence resolving time in positron system, 
t For detector height = 3 can, depth = 5 an, packing fraction = 100 per 
cent, energy threshold = 100 keV for NalfTl) and CsF, 400 kcV for BGO, 
number in parenthesis is photopeak cflBciency. 

Data from Dcrenzo (1981) and Ter-Pogossian « al. (1981). 

Ref. [3, page 46 ] 
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C3iapter 2 

RECONSTRUCTION METHODS 

In general terms the reconstruction problem is to determine 
the distribution of a function f(r,<}5) from its projections p(t,0). 
These two functions are related by# 

p(t,e) = / f(r^<p) ds (l«l) 

c(t,e) 

where c(t^6) is a straight line that has a perpendicular distance, 
t / from the origin and the normal to which is at an angle, 0, 
with respect to the x-axis, fig. i2)» This problem has, of 
course, been solved and forms the basis of transmission computer 
tomography TCT* 

In PET the unknown f\miction is the activity distribution. 
Because of attenuation of the emitted gamma raysi the projections 
p(t.,©) are not directly measured. Rather, the attenuated, 
projection p_(t,0) are measured* These are related to the 
unattaiuated projection p(t,0) by 

p_,(t,e^ = p(t,e) pj^(t.0) 

where p,,(t,0) = exp [“* / U(r, (p)ds] 

^ c(t.,e) 

)Ji(rf(p) - attenuation co-efficient and c(t,6) is the 


same as above* 
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Source 



FIG. 2. parallel -BEAM GEOMETRY. 
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Hence if is known^ the can be calculated* 

The activity distribution f(r, (p) can be computed using the same 
methods employed in TCT or by any method* 

^ ^ ^ ^o^voi'^'tion Back Pro jection Method : 

The data from the PET scanner is the attenuated projection 
The projection data p(t, e) defined by 

p(t.#0) = / X(r,cp) ds 
c 

is related to p^Ct/S) by 

p(t^6) = Pj^^'t^0)/p (t,0) 

where, 

p^(t,©) = exp [- / ds] 

jtx = absorption co-efficient 
K = emission density 

t = distance between origin and the path of 
radiation (RAy) 

0 = angle between xr-axis and the perpendicular 

on the path of radiation 

c = chord (path of radiation) with centra as the 
point (t cos 0, t sin 0) 

p(R, 0) denotes the 1-D Fourier transform of projection with 
respect to s, at an angle 6, 

$(R,e) = / p(t,0) dt 


(2*3) 
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The projection theorem for Fourier transform or projection 
slice theorem [l] , is 

p(R*e) = jd(R cos e, R sin e) (2.4) 

Taking the inverse Fourier transform of (2*4) we get. 


u.(x,y) = 7 “ p(R,e) cos e + y sin e),^, ^ (2.5, 

o *■«> 

For the computation of the integral involved in above 
ec[uation* we need to make scsne approximation because limits 
-oo to oo on R introduces divergence* Therefore., a function 
w(r) is introduced in eqp* (2*5) called as window function 
or filter* Window function is chosen in such a way that 


w( R) 


some function of R* iRl < R^^ 
O IRI > R^ 


Where r = i As, 

V-* 

As = distance between rays* 

^ftar introducing the window function in equation (2*5) we gat 


txCx,y) = 7 z'" p(R.e)w(R) ®’lRl dR 36 

O -R„ 

( 2 * 6 ) 

Substituting for p(R*©) from eqn* (2*3) and interchanging the 
order of integration over s and R we get 
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n «> 

« / / p(s,e) q(x cos G + y sin 0 -- s) ds d© (2.»7) 

O ••oo 


If we assiame object to be of xinit radius then limits on s will 
be frcxn -i to +1, so ws get 


n 1 

M(x/y) = / / p(s,G) q(x cos 0 t y sin 0 ~ s)ds d© 

0 “1 

where, 1/2 As 

^ ^ ^ ^ ® ^ O TrD o 

q(s) = / iRlw(R) dR 

-1/2 As 


( 2 • 3 ) 


(2.9) 


q(s) is called the convolving function* 

“Ihe equation (2»3) and (2# 9) form the basis of the method 
of convolution back projection [?]• If we break eqn* (2*3) 
into two equations : 


p(s',0) = J p(s,©5 q(s'-s)ds 

-1 


u,(x, y) = 1 p(x cos 0 y sin 6, 0)d0 
o 


( 2 . 10 ) 

(2.11) 


The operation represented by (2*10) is convolution and 
the operation represented by (2.11) is knovm as back projection. 
Thus the convolution back projection method applied to PET 
data consist of the following st^^s 1 

( 1 ) Estimation of p,,(t,0) by using transmission tomography 

( 2 ) Obtaining p(t,0) by the relation 

p(t,0) = P-^(t^0)/p^(t,0) 
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(3) Application of convolution back projection method to 
p(t / 0)-. 

E M ALGORITHM FOR PET : 

EM algorithm is derived [ 2 ] frcm the fact that each 
iteration of the method consists of an expectation step followed 
by a maximization step. The physical and statistical considerations 
that lead to a mathematical model of PET are discussed in detail 
in [4] ^ [6] and [9]^ 

The EM algorithm for PET has much in common with other 
iterative algorithms for image reconstruction from projections* 

It is an example of the finite saries expansion reconstruction 
methods reviewed in [l] and discussion of its place, in the 
overall scheme, is given in f5]* We make the usual assurtption 
that the image is subdivided into many small picture elan=nts 
(pixel si-* 

Let X(bl be the (unknown) number of positron annihilation 
events that occurad in the b pixel (1 for which the 

corresponding pairs of photons were detected by the data 
acqiaisition systam* The nxamber of photon pairs detected in 
coincidence by the d^ pair of detector elements is denoted 

by n*(d) (1 < d < D)* The probability that an event in the 

th t-jhi 

b pixel is detected by the d detector pair is denoted by 

p(b,d) while p(b) denotes the probability that an event in the 
til 

b pixel is detected by the data acquisition system as a whole* 
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Ther'Sfore, 

D 

p(b) = E pCb,d) , d = (2.12) 

d=l 

Since the gamma rays are emirted at randOT angles in tlire 
diman S'. ons/ no practical detector syst^^m can intercept all of 
them*, H-i-nce the value of p(b) is less than one for all b* 

Now consider further the time coincident gamma rays 
detected by the system as a whole for which the corresponding 
vevents occur ad in the b pixel* Of these, 1st p(b,d) denote 
the fraction that is detected by the d^ pair of elements in -th 
system* 

By definition, 

. p(b,d) . , 

p'(b,d) = (2*13) 

D 

and E p'(b,d) = 1 (2.14) 

d=l 

Using the following notation for the inner product. 


<p(d), X> = E p(b,d) X(b) (2.15) 

b=l 

the EM algorithm derived from the PET model [4] , [^6 ] and [9] 
may be written as 


= X(b)^ S 
d=l 


, n*(d) , 

( ) 

<p(d),X > 


p'(b,d) 


X (b) 


(2*16) 
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At each iteration (indexed by k, k > O) the algorithm 
produces a new image vector from the old image vector 

The initial image vector k^ is chosen to satisfy^ 

k^*^\b) > 0 , b = 

Due to multiplicative nature of the iterative step any pixel 
of zero value will remain zero in subsequent iterations# 

Two further properties of the iterative method are worth 
noting because they are particularly useful in PET : 

(1) Given a non-negative initial image non— negative 

p(b,d) and n'^(d) (as is the case of PET), all images produced 
by the algorithm are non-negative# 

(k) 

(2) For each image produced by the algorithm (i*e» k for 

k > l), the sum of events in the image is equial to the sun of the 
counts in the data# In other words, the total number of E k(b) of 
estimated counts is equal to the total number Sn*(d) of 
observed counts at each step k , k > 1. 

The algorithm is therefore seljf- normal! zing with each 
iteration producing a redistribution of the activity in the 
image, without any net increase or decrease in activity# 

From the practical point of view, this property and the 
non- negativity of the image may be more inportant than the 
fact that the algorithm proceeds in the direction of maximizing 
the likelihood function* 
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The main disadvantage of this algorithm is the slow rata 
at which it iterates toward the acceptable image* But there 
are soma modified version of EM algorithm which changes the 
image by a larger amount at each iterative step and also retains 
the non-negativity and s el fr* normalization properties of the 
original algorithm* 

CHOICE OF p(b*d) : 

The EM iterative procedure (2*16) estimates from 

old 

an old estimate X by 


new 


(b) = 




D 

Z 

d=l 


n'^(d) p'(b,d) 

b 

2 p(b'.d) 

b'=i 


i ^ ^ B 

(2*17) 


Since we are interested only in those photons which have 
been detected in any of the detector tube^, the conditional 
probability p'(b^d)^ can be taken to be equal to p(b/d)* 

According to Shepp and Vardi's assumption [8^* a positron 
til. 

emitted in b pixel immediately runs to the centre of b and 
then travels a distance p with probability density 

p(p) = pAr/cr^-P^)^ , 0 < p < R, (2*13) 

in a random direction, to reach a point Q in the disk of radius 
R, about the centre of b* Upon reaching Q the positron 
annihilates and the photon pair chooses a line through Q, not 
at random, but uniformly over the n tube direction* 
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Based on the above assumption p(b/d) can be calculated ss/ 
p(b,d) = (2nR)“^ * (width of the intersection of the circle 
of radius R about the centre of b and the strip defined by ths 
tube d), 

where R = radius of the inscribed circle to b 

n = no* of detector elstients around the detector circle* 


Simulation : 

For simulation Shepp and Vardi [9] use the phantom with 
density X(x,y)* 10*^ points are then chosen independently frcm 

the density x(x,y)///X 4 Each point is chosen by talcing a 
uniform random point in the square 1x1/ lyl < then 

accepted with probability X* After choosing each point they 
choose a random line L through the point and increment the 
count for the detector unit corresponding to the two detector 
intervals that L passes through. This gives the numbers n*(d) 


d = 1 f 


MAXIMUM LIKELIHOOD FUNCTION OF X GIVEN n* : 

The variable/ nCb/d) the number of emissions in box b 
detected in tube d* are independent Poisson variables with meen 

E n(b/d^ = X(b/d) = X(b) p(b/d) (2.19) 


Therefore the lilcelihood fxjnction is given by. 


L(X) = p(n*l X) 


E IT g,*"X(b/d) X(b/d)^^^*^l 

A b=l...B n(b/d)! 

d=l ♦ ♦ » D 


( 2 . 20 ) 
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where the sum is over all arrays A of n(b,d)'s with 

def B 

n*(d) = n(»,d) == S n(b,d)^ d =s (2*21) 

b=l 

”tiln 

observed counts in the d detector tube* Also 
def ^ 

n(b) = n(b,,) == S n(b,d) (2*22) 

d=i 

*t-]n 

is the { ■unobservable) number of emissions in the b box» 

According to Shepp and Vardi fs] the concave natxare of 


1(X) = log [L(X)] 


( 2 • 2 3 ) 


can be proved as follows j 

Differentiating (2- 20) w* r*t* to X(bQ#d) we get. 


3L(X) 

3X(b^,d) 


S { TT e‘ 

A t)—l • ♦ * B 
(irsi ♦ D 


n(b,d) ! 


} 


D nCb *d) 

* E £-1 t 2 j 

d=l X(b„,d) 

Using (2tl9), 


3(L(X) 3L(X) 

* p(b , d)* 

3X(bQ) aX(bQ,d) 


(2.24) 


Therefore from (2.2 4) we get. 
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3L(X) 

= E £ TT 

ax(b^) 

A b=l 
d=l 




-X(b,d) 




n(b,d) I 


D 

S 

d=l 


n(b . d) 

2 * p(b .d)} 


(^• 25 ) 


On assumina that 


D 

E p(b -d) = 1 
d=i 


Sqn, (2*25) 

can be written 

as 



aL(x) 
= S 

TT e“ 

(b,d) 

} 


ax(b_) A 

b=l,* . . ,B 


n(b,d) ! 


sJ 

c3.==i ^ » * f ^ I) 






D 




' ♦ Stb^T 

E 

d=l 

n(bQ,d)l 

(2^26) 


By definition of A on (2*20), the conditional expectation is 
given by 


1 -Kb 

E [n(b )ln\X] = 3t 2 TT e n(b^) 

p(n-^IX)A b=l,..#,B n(b,d)I 

/ « * » / D 

(2*27) 


Comparing the second term in (2»26) with (2*27) and using (2*22 5 
wa get 
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1 aL(X) 1 


p(n*lX) J4. -, 

{-p(n^lX) + — — " E [n(bQ)ln^, X] 


X(bo) 


D 


. . I S E rn(b^,d)l tT, X] 

XTb^T d=i " ° 


(i»2S) 


Since n(b,d)‘'s ere all mutually independent, the d 
term in the sum in (2.28^ is given by 


th 


E [n(b-,d) In*, X] = E [n(b ,d) ln*(d), X] 


n*(d) X(bQ,d) 
- ^ 

2 X(b'd) 
b'*l 


4“V\ 

Putting the d term in (2*2 8) we get. 


(2*29) 


3l(X) 

tf k i 

D 

. V 

n«-(d) X(bQ,d) 

3X(b^) 


' Zj 

d=l 

B 

2 X(b%d) 

b'=l 

ai(X) 


D 

. V 

n*(d) X(bQ) p(bQ,d) 

axTb^ 

^ xl^ 

' 2j 

d=l 

2 X(b') p(b',d) 

b'=i 

3l(X) 

3X(bQ) 

D 

= 1 + 2 
d=l 

n*(d) p(bQ,d) 

2 X(b^)p(b'd) 

b'=l 


(2,30) 


Fran (2*30) it follows that 


1 (\) 

aXCb^) 3\(b^) 


D pCb^^d) p(b'/d) 

d=l B 

[ E X(bM p(bS<a)]* 
b'=l 


For any Z( * % i zCb), (2»31) follows. 


B B 3*^ 1(X) 

S S Z(b_) Z(b. ) 

b =1 b^=l "" aX(b ) 3X(b^) 


(2* 31 ) 


- E n*(d) c^(d^ 
d=l 

B 

S Z(b) p(b,d) 

where c(d) = (2-32) 

S X(b') p(b'd) 
b'=l 

2 

and since n'*'(d) c (d) > 0, the quadratic form is negative semi- 
definite so that iCX) is concave-. 

It has been proved in [2, Th. l] that in each step the 
new estima-ca x’^^^ is an improvement over the old estimate X^^"^ : 

L(X^®^) > (2,33) 

with equality if and only if 

L(X®^^) * maxj^ L(X). 

Therefore, it follows from the concavity of log L(x) that 
all maxima of lCX) are global maxima and so the iteration scheme 
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gives a sequence X# X /•»»#X ••• which converges to a 

global maximum likelihood estimator X • If the maximum of 
L(X) is not unique then by the concavity of log L{ X5 any convex 

combination of maxima is a maximum and X depends on the initiol 
o 

choice X 
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Chapter 3 

Nimerical Simulation and Con elusion 


In our simulation we consider the case when there is no 
absorption of the photonst The measured data are n*(l ) » ,n''*‘( D) 
where n*(d) is the total number of coincidences counted in the 
d detector tube* The simulation of n^(d) is as follows : 

At equidistant points in the detector tube d with \aniform 
spacing 0*01, the positron activity proportional to the intensity 
of the picture, in the pixel, in which the point lies contributes 
to the tube d, data n^(d)*. This is done for 123 xl28 tubes 
corresponding to 123 detectors mounted on the detector ring* 

123 of these cases correspond to a null tube with zero data* 

H 0W 

The iterative procedure for obtaining a new estimate \ 

oX(3. 

from an old estimate a , 






D n^Cd) p(b,d) 

£ — 

d=i 2 
b'=l 


is simulated as follows : 


b = i,». *,8 


(3.1) 


Let Tc(b,d) denotes the number of times pixel b has 
contributed to tube d and N(b) is the total contribution of 
pixel b to the detector ring, then p(b, d) is given by 
k(b,d) k(b,d) 


p(b,d) = 


S k(b,d') N(b) 

d' 


(3.2) 
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and wa have 


S p(b/d) = i (3.3) 

d 

Now if Pic(b) is the value of pixel b^ then 


X(b) « PlC(b) * N(b) (3.4) 

and 

k(b' /d) 

s x_(b') p(b%d) « s Pic_(b')N(b') * — r-rr- 

b' ° b^ ° n(b^) 


= S PlC^(b') k(b',d) 
b' ° 

= nQ( d) ( 3* 5) 

Thar afore (4»l) can be rewritten as 


N(b) PIC;j,(b) 


N(b) PlC^(b) 

NXb-^ 


S 

d 


n*(d) lc(b,d) 


PlC^(b) 


PIC (b) 
Nib’) 


n*(d) lc(b,d) 

E 

d 

o 


(3.6) 


Equation (3. 6) has been implemented with various choice of 
initial approximation* L— 1 and L-2 errors are then calculated 
as follows : 


1^1 ~ Npfec ^ ^^b ** ^b^ * ^ “ I/*' *#6 

b 
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wh ere, 

NPIX = Total no. of pixels considered 
= Reconstructed pixel b 
= Original pixel b- 

Four esses which we have considered are as follows s 

(1) Initial approximation is constant'. 

( 2 ) Initial approximation is from CBP having rebinning 
projection data from PET [PBTCAT^ 

(3) Initial approximation is from CBP [CAT*]. 

( 4 ^ Initial approximation is the last iterate of (!)■» 

Observations : 

(1) Table [3] and Table [ 4 ] shows the corresponding L-1 and 
L“2 errors of each iteration for all the four cases 
respectively. 

( 2 ) Figures of each iterations are in the appendix B* 

C ONCLUSION : 

From the figs* [3-6] of error and iteration number , n, 
we find that as n increases, the curve on the logarithnic scale 
approaches straight line with a negative gradient* 

Hence for large n, the error values on the curve are of 
the form 

}cq^ , where k is a constant and q is a positive number 
less than unity, implying that the converg^ce is of first order* 
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At the outset/ from the reconstructions reported in 
Appendix B, it is clear that the CBP reconstruction in the 
parallel beam set up is generally the best and remains essentially 
unimproved by E-M iterates* The CBP algorithm applied on the 
rebinned PET data itself gives a reasonably good reconstruction 
which is only moderately improved by successive E— M iterates* 
Further starting with a constant initial approximation in the 
E-M algorithm seams to be quite time consiiming* 

The behaviour of errors in the three cases reveals the 
following information about EM-itarations s 

(i) If the initial approximation is crude the subsequent 
errors fall down rather qxiickly/ till we get a good 
approximation of the cross- section* 

(ii) Starting with a good initial approximation the successive 
iterates improve the approximation till we get a close 
enough approximation of the cross- section* However/ in 
this case the decrease in the error is rather slow* 

(iii) If the initial approximation is close enough the following 

iterates do not necessarily improve the approximation and 
the errors, inf act, might even increase* However, this 
increase is rather marginal and the iterates continue to 
remain close enough to the desired cross-section* As the 
pattern of errors is more or less random apparently in i 

this case the errors are solely due to round offs* 
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Combining (i^, (ii) and (iii), it follows that one may 
start rhe procedure with an arbitrary initial approximation and 
continue the process till successive iterates keep on changing 
substantially. One rhen should stop the iterations when the 
changes are only marginal. 

In conclusion^ a practical strategy using EM algorithm 
would be ; 

( a ) Rebin the data 

(b) Apply CBP to get an initial approximation 

(c) Apply Em-algorithm with the initial approximation given 
by (b). 

We may note at this point that rebining the PET data 
for adaptation to parallel beam geometry introduces interpolation 
errors to which it is wall known that the CBP method is quite 
susceptible* On the other hand since the PET data is 
automatically a fan-beam geometry data, the corresponding fan- 
beam CBP may be used. The third alternative is to use iterative 
algebraic reconstruction techniqtie Art, for the PET data. In 
view of the reported superior performance of iterative art 
[ 4, il. 5 ], this possibly should be carefully examined since 
the time taken in EM iterations and Art iterations would be 
comparable, while the non-negativity and constant total activity 
by products of EM algorithm could be incorporated with equal 
ease in the ART case as well* 
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TABLE 3 

; ERROR L-i 




ITER* 

NO. 

15 CONSTANT 

constant 

CAT 

PETCAT 

1 

1,706590 

5, 446335 

0. 31 4174 

1.115359 

U 

1.627746 

4,390611 

0*313572 

1.111345 

3 

1.573573 

3.7232 53 

0. 333133 

1*107734 

4 

1. 517605 

3.274752 

0. 343665 

1*093405 

5 

1.459325 

2.943726 

0, 352091 

1.0312 52 

6 

1.413131 

2*692951 

0. 350537 

1*053933 

7 

1*364129 

2*502359 

0. 349 082 

1*037617 

3 

1. 32 5303 

2*343966 

0, 344 563 

1.019561 

9 

1.2^1002 

2*210653 

0, 342462 

1*009630 

IQ 

1.257399 

2*096299 

0. 34035 5 

0.992981 

11 

1.229010 

1,999097 

0. 341B60 

0.930740 

12 

1.201023 

1*915438 

0. 341360 

0.969004 

13 

1.171833 

1,336394 

0. 343966 

0.9 56365 

14 

1.145652 

1.772796 

0. 343966 

0.943425 

15 

1.120975 

1.714716 

0. 34 5170 

0.932 591 
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TABLE 4 : ERROR L-2 


ITER. 

NO. 

15 CONSTANT 

CONSTANT 

CAT 

P ETC AT 

4 

X 

2,973299 

5.373633 

0.690418 

2,239094 

1 

2» 3-i<±712 

5,777310 

0*651357 

2,133733 

3 

2*770215 

5.1‘±14S3 

0.677433 

2.172443 

4 

2.691199 

4,723436 

0*692160 

2*153457 

5 

2.619115 

4,409243 

0.71272 3 

2.119440 

6 

2.551191 

4,159136 

0.717143 

2.031762 

7 

2*437102 

3.953651 

0* 72 73 51 

2.042433 

3 

2*433466 

3,775326 

0*733234 

2*009307 

9 

2 . 335330 

3.622247 

0*752966 

1.983607 

10 

2.336320 

3.475290 

0*765345 

1.957195 

11 

2.2 51353 

3.359 537 

0.739641 

1,932533 

12 

2.251353 

3.2 54613 

0*309959 

1.910265 

13 

2.207143 

3.153392 

0.334120 

1,376150 

14 

2*169315 

3.067551 

0*353909 

1.354376 

15 

2*129921 

2 9 3Sa» 9i 

0.873574 

1.836124 




Fig. 3. VARIATION OF L-1 ERROR 



Iteration number 


n Pe-t c cvt 
A con s'la n't 
0 Pe-t ca“t 
+ c a*t 


\ 




Fig. 4. VARIATION OF L-2 ERROR 
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TIMS COINCIDENCE TECHNIQUE 


ThivS technique identifies the time rela-fcionship between 
the radiations emitted by a radio-active source* In its simplest 
form^ a doxable coincidence circuit/ when fed with two input 
signals provides an output signal only if both of them arrive 
within a certain time t of each other* This time t / is known 
as the resolving time of the coincidence circuit* 

In an overlap type coincidence circriit accepting two input 
signals each of width T / the resolving time is 2't because 
coincidence output results only if one pulse lies within + i of 
other* Coincidence unit with a resolving time in nanosecond 
(nsec) range is termed as a fast coincidence unit* On the other 
hand/ a coincidence unit with a resolving time in micro-second 
(■]□, sec) range is called as a slow coincidence tjnit* 

By using appropriate detector/ coincidence measurements 
can be made between alpha- and gamma-rays, beta and gamma rays, 
gamma- and gamma rays, neutron and gamma-rays or between any 
such combination* 

Normally the detectors are kept back to back* But by 
keeping one detector in a fixed position and moving the other 
in different angles, we can do angular (coincidence) correlation 
measurements* Such measurements are usefxil in determining the 
spins and parities of excited states of the nuclei* 



The principal problem in the coincidence measuranents are 
the accidental coincidences* To understand this# let us assxjme 
that angular distributions are isotropic# then# no* of 
coincidences detected per sac N # is given by 

O2 

where 

Nq = No# of positron emissions occur/sec# 

= Solid angle subtended by detector 1 
Qj = Solid angle subtended by detector 2« 

If all detected photons are due to these N^^ positron 
emissions# the total count rates in the two detectors and N^ 
are, 

“ ^o ^2 = ^o ^2 

If photons are detected in the two detectors within the resolving 
time of the coincidence analyzer t# they will be recorded as a 
coincidence# The chance that this will happen accidentally 
before or after each photon is received in the first detector 
is 2 t times the probability per second for a photon to reach 
the second detectors# or 2 TN 2 * Since this situation occtirs 
N^ times per second# the total number of accidental coincidences 
per second, # is 

Ng = 2t: * n^N2 . 

Therefore, the ratio of true to accidental coincidences 

is# 




2Mo H 


1 

2 ^ 

o 




For good coincidence measurement/ this ratio should be 
sufficiently large* The ratio can be increased by reducing 
Nq or T • Since reducing means reducing i*e« rate of 
data accumulation it is advisable to keep t as small as 
possible# 



COUNTS DUE TO 
DETECTOR 1 


COINCIDENCE 

COUNT 


COUNTS DUE TO 
DETECTOR 2 
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z:s)tfi3e3Exa2: ibx«:k«k jcatat^rxo 
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raxraxxxxwaxxx icxn w n n «i-ica ^cxxxKKXMXicxxifEtsiax m -4 

HI>ihOOOOO XCdDXXXXXacXXWXXCBZXOOM -I- 
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aaaaaaaaa£^^MM 44 4 MJ 3 3 /)^^a aM» 4 44 4 4 444 M.oaaMa aaa- 44 a a as 
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xcox—iaBaaaaaxaxo MM- az oxxxzcom* M>i >4444 4 «‘Co Mcotfsaoxaaaaaaaa 
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CB SCH B B ttK » as & n (C 2 H- u at K Bib KB VB BBSat XK S v: 

BBXBBBgBBBBKB^W i:aiBB)K-M>«> -(-BBBBBBBBBiCBBSEK*; 

>((BSBBBBBBttBBB:e:^ #2 •-4' H aE««r)«BC: B 2 MS ? 

BSarBBBjBBBBBBBadD^- | «03t-»»- H II ««*■«*«■» IX oa* V3 

■BBBBBBBBBBBBBB X BB3K2 'i>2 BBB ZBBBBBBBIB S 

BBBBSSBBffiflrBBaBBBBSl mvi% il K> COIBB^K il Ml SCaCBBIBZZ J3:S *^2 M 

*;S255E!ISSIS25*22*** . -» + « « urBiBBBBBKBKa^irBatfcizcGC mn-* 
aaBBaEBBBBBBBBBBBBUI >31 J3BB»B»«B£XBBSBBSBSXaB33I OtOil Halt 
BBBBBBB^BBBBBBBBBB^^^^-BBt/^BBBBBBH-CH-S^aEBScaRBBBaBBIBSlCtl tCBt iC 4- 
KJEBBBBBBBBBBaVBBBBlKfilBIBBaaMIS; Cl» 3E®IW2aa«OH ZW*«M<s<O02u-,|f u « 
.BaCBBBBBBBBBBBBBBBKBBBEICBBBBBOCt-lCHZZ It ■c;:3rB23»K»3 xnv.^l BSIC 

3CtE?BBB3BBOaE:BaBaaS!SBBaiBB&BraaB»2»tlBattBflS3K HI tHlI »-(»•)(■ t-m-4H4 0XZ0HigB}4(j3K;tf3)-<4 II «f 

BBBSBBBB BBBBBBBBBBBBBBBBBBBBBBBBBBBtBBKtO !/:«-»- II V3C/: 3c 2 2 2 tfcKKIZ cCt/3? 
M<aCBBSSBaBS3i&iBaiSBBaiBBBBBSBBBBBBBBI»»giBBaBBXaa4>4‘«- 4-»:SZ3--3tzaz2V2Z il a 
.BBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBIBKKBBBCB* 11 4^ 4^ 4- 4'U: U 21 tKKX 21 MC 
HBBBBffllBBBBBBBSBBSBaBBUIBIBBSrBBB»BKIX3eiZ:2ZZZ3B4>4>4-4-4>U3u:iua2XHaaMC0&0 
SffiBBBBBBBBBBSBBBBBBBBBBBBBBBBCnc/^C M «- » «> 4*- 4> 4- 4- «■ Xr«- 4- 4- 4- 4-)- H i-^ZtES 3C2BXIIB 
SeSBaBBflBBaSBBCVBBBBBBBBBBBaaBaBSBtimt U tt«' 4>»*-44‘4-4^*«03Ht4-4>4>44>^O3;H«aiZa® 
>»^BBBBBBBBBBBBBBB»BBBBBBBBBBB»-<23 It »«.4-»««- •• II H l/lMOl 4- » 4^ -f 4* 4-» CX S t£ 2 S BB IK 
SaSiWBBBSDBBBlimBBBaBaiBBBBBBBSBZaEIlKOMll «- » 4>» » 4- 11 H 11 H Kit H il 4> 4- 4> 4>» J3'J3Z I KB XBB 
BBBBBBBBBBBBBBBBBBBBBBBBBBB2I n |tcCH<t/:}«H-4i«- 4'4> It U II H Z II 4' •• tl 4> - li ih* tCB KBXB 
KBaBBBBBWaSBBBBBffiBBBaBSBBBBBZ: R M>^H i-(ll 4-4-4-11 il 4'X4-I1 4- ••It » I ••«-«- r4Z3B BK KB BK 

IBG[ aBBBBBlBaBBBBBBBBBBiBaBBBBBBBB3»-^H4K 4'4(> 4 4- •« M 11 li 4 it It »1 It •« |l 4 4- 44 M 2 2 XB IBiBB 

COXBBBBaSSBGBBBBfiSBBBilB0aB»aaBSBZ4 44 4 44 44M44 4 11 II il i| || II 11 4 II ••» 4 4 44 4 O H Z X mBBB 

CdLBBBBBBBBBaBBBBBaaBBBBBaB2)>4 4 •• II 44 l~ 44 4 4 11 M U 11 il ti U II tl U •• It 4 4 4 4 4 C^4 C X BBBB 

nsBaBBBaaBaBaaaBBBaaaaaaaBB>^4 4 4**>*4v^i'^H(4 n 4 ti it it » it ti it iiti ii ti it 44411 ii 4 ti nzasaa 

CXBBBBBaaBSBBBBBSBBBBaBBBaBCt»4 tt •• I II 4 4l-F^4 Ii 4 4 II It •• il 11 ft It •• 11 It 4 4 4 It ••*.•«•• 40XmKB 

z®aaaaaaara:araassassaaaBaaaaazc04 u 4>«zst 11 ■aco44 4444 4444*^(44444411 4 ii 11 ii aiKaa 

C»:BBBBBBBB4K.2 2«rBBBB;»BBBBBB4 » I •« 4 I |i 4 4 4 •• It |{ |{ |i li 1 f I 8 I 4CSIEB 

2'saaa!KasssrairasssffBBaaaaaaaaB4 it •~t 4 t/3i-(4 4411 m it 11 11 ti 11 it ii 4 it 444 44 ii •••••• -4 hbb 

ii -4 4 4 1-^4 1/2 v:tC4 4 4 II » It U - 4 II ••4 444 R 1-^4 4 •••42tRBB 

SXSEffl«a28ffiBaSRrBaSBaffia«l3SSa«-il»‘C442CC3CiiK4»44IIUtt»ll44M4CttH»4-*.|4ZaiS« 

K 11 ••4 tl tf MH-tr.BZtCOK 44 411 4 tl 4 4 4 4H4|t ft ••••-4<€B«: 

X X KBfflBaBaaBaxaaatn a n4ii{icti ii^ 43tta2accu4 4tt 444444(-X444aa = »- **4a8 

{T IS:BiBBBBS^BaBBB»BBaaBBBaBa:4 4 11 B If 11 4BB2tC2ZXX » ^4 4 H 4 4 4 »h>^U BK4 4 C? Z 4 4 • ' ' 

a = »aafflra aBaararvaaaaaaaaa ii4it««K42xactiic^'-/’22'./’B44444»-^4Raa>-4‘^ax>-c/44aa 
B(/:i«.SILjS!BBBBBBBBBBBBBBBBC 4>-t4 4B2 It BBZ ttKZ 2ZZtC»H>.< 4 4 4 44 ^h-«4BBB h«CB CX Z 2BBB 
>*[CC«aaeaaaaaeaaaB8ax3ti4<f«4ax4aaa44aa2sirtf444444H-44ac«ac:23z3asx« 
Zt/:S:S:SlBaBiBaBBBBBBB2 11 U II M4 U 4BX:2BBffl:4eiH IBUU 44 444 4 4 4 4 4 aBXBBC^ OBBWBa H 
•!a:NxaaaaaaaaaaaBx««««-N4ii«sa2a8aBaifa2Ct-44H32a444»-aaaxas2asaaa 

]4t/:C>KBBaBBBiSBBE3iaBa'4 El ••tEBZlBaiE»cCBBCiC4 44 U B 44 » 4 44BBIiB:B3 fEBBBBIK 

2ii3aaaaaaaBaaaa2-u«**<-*<i laBasaa^iBasriiif fr4tiiiKii4«n4B»kaaaa49aa«B 
XtOt/^BBBBBBBBBBBKZ R 11 •*••11 it •• ••XBfiCV^IBBBV:: II 44 4 tl 4 4l-t4)K K- 11 H-BXBBICKnCBBSBB 
4{r2isaiBaaffiBSBa3irM»-*«ii44«^^ 4aaMiH3M244 4i(»i ti••u4MB rtU3a2a*aBSBisaB 
4CV:XBBBBg3SBBn»XCC: il li II 4 411 •• lIBaZ Bt/:N 4 U n •• tl 11 R It II tlBtl RaCBZlEKBBBBBBBB: 

M>cr2aaaeasaBa3ac/ HU(-ii4 44 ti i4iaBax2i->i 11 11 11 tl••44lllt^ HaasaBsaassiaB 
SEbCBBBBDBaaarZ znuu 11 44441 ua:BBB{S:i>~«*ll 11 11 it tl 11 4 tl U 11 t/ 2 B 2 B:aCBBKBBBB 

cOZ:ZaSSBB!BaB!SB:^v3ll 114 114 4 4 4tl>HB^MB4 II It tilt It II 11 It It tl •• fl BaXaBBaDBaBSSB 
,CZaB.BBBBBBB:a;Z>-bC tl ll»H4 4 44«>4tCtf:i24 4KIC)KW«* 114 il II It ILBXBBBKBiCBB 
accBaaaaaaBa!/3zii It 384 444 4tojacoj04BaaBaBtfiz4 •! 1 aeatBaBaBBsa 
Ctt:cBaaBBBaBit/:ci:^R t^(-<4 44H<c^Bict/:ts:arBaaB>K2CKBcA4BXBBBBBaa 

zB3Eaaaia3aB3coao3t/3» a 4N4co33cs:i£aaaazaacozB3;BaBaBaBBaa 
BSIBtSBBBaiBZZeCtca a U RO^CCZ U ZBBBBZBCC CCC &BBBBSEB»aB 
a a as 32 BBa B K to a a o '/la B z ZB 22 ma 96 3 : B a u3 33: BB BB BBa I 

a a a B a B B X 2 3 c to 2 (t 2 a GT X z «: >K 2 c KO B tx CO c c c RF B » K B a 
3D at be a a a 33 a B K K I >« Mt X B >« to a z o a 'J3 a z z o a to I a B a 2 
XB B B K a B B a s: X IB ss X c 2 z dc cooo o CO c CIO asaas 
BaaBBaBas3aB'J3'J3oz2ioooaz z2aaa; 
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DO 232 Ksl,i2 
PHOJ(t‘0=SOHT(.F{iOATC32-K)) 
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CONTINUE 
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DO 901 J=i;.'IPIC 

IFCPICNEW(i,J3.NE»0.03 GO TO 9 
PICMEWCI, J3SCODEC32) 

GO TO 901 

K=C32.5-PlcaLDCI,J)*PlCOLD(I,J)) 

XFCK,1,T,13K=1 



f 


^i}i 


543 

SOi 

576 


IF(r,aT.3/.')h=32 

PlC'.f!:;>''Cl,d)=CriDF(K) 

r 3 r f ! p 

on 543 1 = 3 . IPIC 

.VRJT!!:C24,>) (P'^C^£.'^( l, J) , J: 

CO..Tlf;L'.“; 

C:1 


V. 

STHt 

f::'D 


TTi'rJF 






